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In this paper we present an exact Grassmann stochastic Schrodinger equation for the dynamics 
of an open fermionic quantum system coupled to a reservoir consisting of a finite or infinite number 
of fermions. We use this stochastic approach to derive the exact master equation for a fermionic 
system strongly coupled to electronic reservoirs. The generality and applicability of this Grassmann 
stochastic approach is justified and exemplified by several quantum open system problems concerning 
quantum decoherence and quantum transport for both vacuum and finite-temperature fermionic 
reservoirs. We show that the quantum coherence property of the quantum dot system can be 
profoundly modified by the environment memory. 
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Introduction. — Quantum dynamics of electronic sys- 
tems coupled to fermion (electrons) or boson (phonons) 
environments has recently attracted wide-spread inter- 
est in quantum transport, quantum computing, and 
nanoscience [l|, y]. For example, the size reduction of 
quantum devices in microelectronics requires controllable 
systems consisting of only a few electrons, where quan- 
tum coherence and quantum interference become dom- 
inant. In addition, quantum dots coupled to electrons 
of a metal is an interesting setup in quantum informa- 
tion processing where the quantum coherence of c^ubits 
is essential for many quantum algorithms to gain an ex- 
ponential speed-up over their classical counterparts. In 
the case of quantum open systems immersed in a bosonic 
environment, a theoretical formalism for quantum open 
system dynamics was developed to provide a powerful 
tool in studying quantum systems in a non-Markovian 
regime [3|-l9|. This theory has extended the standard 
Markov quantum state diffusion [lO| to the case involving 
strong system-environment interaction or structured en- 
vironments. For a Markov fermionic environment, both 
the quantum state diffusion equation [10| and Lindblad 
master equations llj can be used to describe quantum 
dynamics of the system of interest. However, for a generic 
non-Markovian fermionic environment, a stochastic the- 
ory analogous to the non-Markovian quantum state dif- 
fusion equation is still missing. The major theme of this 
Letter is to develop a non-Markovian stochastic theory 
of electronic systems coupled to a fermionic environment. 
As an illustration of the power of the stochastic approach 
developed here, we derive the exact master equations 
governing the reduced density operator of the electronic 
systems coupled to vacuum and finite-temperature reser- 
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To begin, we consider a model involving an electronic 
system in contact with a single fermion reservoir, where 
the system anti-commute with the bath |12l ]. The gen- 
eralization to two reservoirs {e.g., source and drain) can 
be established in a similar way. With necessary modi- 
fications, the formalism is versatile enough to deal with 
stochastic gate potentials and nonlinear couplings. The 
total Hamiltonian for the system plus environment may 
be written as [l[. 



Htot =Hs + Hn + Hi, 



(1) 



where Hg is the Hamiltonian of the electronic system in 
the absence of the environment, H^ is the Hamiltonian 
for a fermionic reservoir; Hb. = Ska ^'^k^ka^ka where 



^ka ' ^ka are the fermionic creation and annihilation oper- 
ators {^ka, &1(;/q,/} = (^ka.k'Q', and the interaction Hamil- 
tonian Hi is given by 



Hi = fl'^itkaLbl^ + tkabwaL''), 



(2) 



where L is the system coupling operator and tka are the 
coupling constants. Note that Hg is an arbitrary Hamil- 
tonian operator that may contain interaction terms for 
the system particles {e.g.,, Coulomb interactions between 
two electrons). The coupling operator L may, in general, 
be represented by a set of fermionic operators which are 
coupled to all the participating external agents such as 
the source and drain reservoirs. 

The primary purpose of the present paper is to de- 
velop a systematic theory for the models described by 
([T]) and ([2]), that are relevant to non-equilibrium statisti- 
cal mechanics. Path-integral theory and quantum devices 
based on quantum dots and mesoscopic electronics |13l - 
ligj l . In the framework of the stochastic Schrodinger equa- 
tion (SSE), the state of the open quantum system is de- 
scribed by a stochastic pure state, which is generated by 



a stochastic process. For the fermionic environment, sim- 
ilar to the fermionic path integral, the fermionic stochas- 
tic theory will involve a Grassmann stochastic process. 
We show that the reduced density matrix of the system 
of interest can be reconstructed from the pure states by 
taking the statistical mean over the Grassmann noise. As 
such, in principle, the exact master equation governing 
the reduced density operator can be recovered from the 
SSE, as illustrated by several physically interesting mod- 
els below. Compared with the bosonic case, it is more 
complex to recover the density matrix from the corre- 
sponding SSE due to fermionic features of the system 



and bath 0,[21| 



Non- commutative stochastic Schrodinger equation 
(SSE) and master equation. — Consider the model de- 
scribed by the total Hamiltonian Eq. ([1]), in the interac- 
tion picture with respect to the fermionic reservoir it be- 
comes (setting h = 1), i?/ot(0 = Hs + iJ2j t,Dhje-'"^^^ + 
h.c), here the subscript ka is suppressed as j. In or- 
der to trace out the environmental variables, we intro- 
duce the fermionic coherent states |<^j) which is defined 
as bj\i^j) = S,j\i,j)- Here, S,j is a Grassmann variable, 
s^isfying {S.^,S,j} = 0, and {S,j,bj} = fe,6]} = 
[22| . As shown below, the derivations and results for 
fermionic SSE are more complex than the bosonic case. 
Using the fermionic coherent state, we define |V't(^*)) = 



(^|e»ffR* 



iHu 



*|^tot(0))where|0 = |6)®l6) 



IS a 



collective notation for the coherent state of the fermionic 
bath. Then, the non-commutative stochastic Schrodinger 
equation for the system of interest can be written as [2l| 

dt\i>t{C)) = -iH.slMn), (3) 

where the effective Hamiltonian is 

iJcff ^Hs- iLCb t-iL^ f dsK{t, s)t ^> (4) 

where ^l^ = —'iJ2k^k^^'^''^^bk i^ the Grassmann noise 

and the correlation function is K{t, s) — J2k as* a£*'° ~ 

J2k |ifepe^*'^'=^*^''^ Note that our approach is applicable 
to arbitrary correlation functions for both Markov and 
non-Mar kovian cases. 

Eq. ([3]) may serve as a fundamental equation for open 
electronic systems coupled to a fermionic environment. 
Crucial to the practical applications of Eq. (|3]) is to 
express the Grassmann functional derivative under the 
memory integral in Q in terms of system operators (^ 
0]. In order to calculate the functional derivative in the 
stochastic Schrodinger equation, we introduce an oper- 
ator called the fermionic Q operator (similar to the O 
operator in bosonic case) as. 



Qit,s,n\Mn) 



;\Mn), 



(5) 



where the functional derivative is alwarys assumed as left 
functional derivative throughout this paper. With this 



definition, the effective Hamiltonian in Eq. Q can be 
written as, 



^cff' = Hs- iUl,,t - iL^Q, 



(6) 



where Q{t,$,*) — /^ dsK{t,s)Q{t,s,^*). The stochastic 
equation is derived directly from the microscopic Hamil- 
tonian without any approximation. It should be empha- 
sized that the system Hamiltonian Hs and the coupling 
operator L are entirely general. The evolution of the 
electronic system is governed by this non-commutative 
stochastic equation. Although the mathematical form of 
the equation is similar to the bosonic case, the behavior 
of the system can be different due to fermionic features 
of the environment [12| . 

From the consistency condition for the fermionic SSE 

5^.dt\ijjt{^*)) = dto^^jipti^*)), the fermionic Q opera- 
tor satisfies the following equation. 



dt 



= -i[H,ff, Q] - «^€* (^eff -Hs). 



(7) 



Once the fermionic Q operator is determined, the SSE 
can be cast into a time-local stochastic equation with 
the Grassmann type noise. Note that the reduced density 
operator for the open fermionic system can be obtained 
by taking the statistical average over all the Grassmann 
trajectories. 



Pr = jVg[i]P 



+ \MC)r{M-£.)r 






(8) 



(9) 



where 2?g[^] — Yik'^^k ' d^kS^^'''^'' is the Grassmann 
Gaussian measure. The notation "~" stands for the 
transformation that change all the creation and annihila- 
tion operators to their additive inverse, and " "'^ " means 
the transformation " ~ " is applied to the system of inter- 
est only. Note that the reconstruction of density opera- 
tor in this sense is more complex than that in the case of 
bosonic bath [20|. Due to the anti-commutative features, 
the density matrix is composed of four different kinds of 
trajectories in Eq. ([9|). However, they can be generated 
by the same trajectory |i/'t(^*)) via transformations [21| . 
so we will only consider \tpt{^*)) in the following parts 
and use the shorthand notation lipt) = |V't (€*)}• Taking 
the time derivative and using the Novikov-type theorem 
for the Grassmann noise [2l|, we can derive the formal 
exact master equation as. 



dtpr 



-l[Hs,Pr]+{JVg[^][QP.L^]+h.C.}, 



(10) 



The above derivation of SSE and master equation is only 
for the case of vacuum , in which we assume the sys- 
tem and environment are initially in the state |V'tot(0)) = 



|V'(0))s ® |vac)R. However, the finite temperature can be 
incorporated in our approach, as shown in the examples 
below [21]. 

In a special case wherein the fermionic Q operator is 
independent of noise, the master equation Eq. ([TU| be- 
comes a simpler form 



dtPr = -i[Hs,Pr\ + [Qp, L"^] + h.c. 



(11) 



If we take the Markovian correlation function K{t^ s) = 
T5{t, s), the Q operator reduces Q = TL/2, and the mas- 
ter equation will reduce to the standard Lindblad Markov 
master equation, dtPr = — *[-ffsjPr] +r/2[ijOr,L^] -I- h.c. 
Many-fermion system coupled to a vacuum fermionic 
reservoir. — In the first example, we consider a many- 
fermion open system coupled to a fermionic bath initially 
in the vacuum state. The total Hamiltonian is i^tot — 
'E'jLi ^jd]dj+J2k ^kblbk+J2jM t:i.kdybk+tj,kbldj, where 
dj and d, [j = 1 to N) are the fermionic annihilation and 

creation operators for the fermions in the system, bk and 
&J, are the fermionic annihilation and creation operators 
for the bath. Here Ho, ~ ^j^jd.dj, and the coupling 
operator ^s L — '^- dj . Then, the fermionic SSE is given 
by 

^dti^t) = {Y,^,d]d,~iY,m,~iY.d]Q)\4>t), (12) 



where the Q operator is given by Q = X)i fii^i s)di. Sub- 
stituting the Q operator into the Eq. ([7]), we obtain the 
differential equation for the time-dependent coefficient 



f,it,s) 



Mt, 



injfj{t,s) 



^N 



rLiFk{t)f,it,s) 

with the initial condition fj{t, s — t) — 1. Fj{t) is defined 
as Fj{t) — /p K{t, s)fj{t, s). Thus, the exact Q operator 
is fully determined. Then one immediately obtains the 
exact master equation 



dtPr = -i[Hs,Pr] 



N 



N 



^h.c.}. (13) 



We have shown that the SSE can be a power tool in 
deriving the exact master equation for an exactly solvable 
many-body system. 

Single quantum dot (QD) coupled to a finite- 
tem,perature hath. — We consider a single QD interacting 
with a finite-temperature fermionic bath. In the standard 
Hamiltonian Eq. ([T]) and Eq. ©, i?s = oj^d^d and L = d. 
It is known that the finite temperature model can be 
transformed into the vacuum case by introducing a ficti- 
tious bath "c" as [a,|9|, Hq = — J2k ^kCkCk- By properly 
choosing the parameters of bath "c" , the combined bath 
"&-|-c" can be initially prepared in a pure state which can 
be transformed to vacuum, while the real bath " &" is pre- 
pared in the thermal state after tracing over the fictitious 
bath "c". I.e., pb{0) = Ti^pbciO)] = exp[-%^]/Z 




FIG. 1: Time evolution of the coefficients for single quan- 
tum dot in a finite temperature bath with diflterent band- 
width. The real (imaginary) part of the coefficients Fi {F2) 
are plotted separately. The parameters are T — lOOmK, 
/i = 2 X 10"^ey, a;o = 3 X IQ-'^el/. 



where Z = Trcxp[ 'l^tt ] is the partition function. 

Therefore, the finite temperature model is transformed 
into an effective vacuum case with a fictitious bath |2l|. 
Then, the exact SSE for the QD system is given by 

idt\tpt) = {uJad) d ^ id) Qb' - ikl' t + ^d^Cc' t + J^QcOl^t) 

' (14) 

where the fermionic Q operators are defined as 

Qb'li^t) = J^dsKb'{t,s)'ti'^,JiJt), Qc'Ht) = 

J^dsKc'{t,s)6^'^^ J-^t). 

Following the scheme of deriving master equation from 
the general SSE, we obtain the exact master equation of 
this model, 

dtPr = -i[Hs,Pr] + {Flit)[dpr,d^]+F2it)[prd,d^+h.C.}, 

(15) 
The dynamic evolution is determined by the time- 
dependent coefficients. In the Markov limit, the coef- 
ficients are time-independent constants, however in the 
general non-Markovian case, the coefficients will be af- 
fected by the parameters of the environment. In Fig. [T] 
, we plot the dynamic evolution of the coefficients Fi (t) 
and F2{t). For simplicity, we choose noise-free Q opera- 
tor in our numerical simulations. The spectral density is 
chosen as the Lorentzian form i^.(a;fe)Aa; = — — ^h -,2_i_i,2 ■ 

When the bandwidth b is large, which corresponds to a 
white noise, then the coefficients conv erge to constants 
rapidly, approaching the Markov limit [l8| . On the con- 
trary, if the bandwidth b is narrow, the distribution of 
the spectral density represents a colored noise, then we 
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FIG. 2: Time evolution of pn with different parameters of the 
bath, (a) is plotted with different temperature T, and (b) is 
plotted with different bandwidth of the spectral density. The 
other parameters are fi — 2 x lO^^el/, ojq = 3 x lO'^el/. 
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expect that the non-Markovian properties (e.g. time- 
dependent coefficients) will be observed. As a result, we 
see that in Fig. [5] (b), the density matrix element pn 
converge to the steady state faster when the bandwidth 
b is larger, and it exhibits significant fiuctuations when b 
is small. 

Another parameter that will affect the non-Markovian 
properties is the temperature of the bath. As it is shown 
in Fig. [2] (a) , the non-Markovian behaviors become ap- 
parent in the low temperature regimes and the Markov 
dynamics at high temperature. 

Double QDs coupled to two finite-temperature 
fermionic baths - In the third example, we consider 
an electronic system coupled to two fermionic baths 
described by the following total Hamiltonian, 

Htot = ujid\di + U!2dld2 + gd\d2 + g*dldi 

k 

+ {Y^ t2,kdlb2,k + Y^ Kkd\Kk + h.c.}, (16) 

k k 

where di and d\ {i = 1,2) are the fermionic annihilation 
and creation operators for the two quantum dots in the 
system, and bi^k and &| ^ are the fermionic operators for 
the baths. This Hamiltonian describes a physical model 
that double quantum dots coupled to two fermionic 
baths, the "source" and the "drain". This model is 
widely studied during the past few years with the input- 
output theory and fermionic path integral [19i, ,23] . 

In our SSE approach, following the similar procedure 
in example 2, the effective Hamiltonian of the model in 
the exact SSE can be established as 

idt\i^t) ^ {Hs - idlQi.b' - idiS,l,b',t + id\^i,c',t + idiQi.c' 

- id\Q2.b' - id2^2,b',t + «4C2,c',t + «^2Q2,c')IV't) 

(17) 
where 0^,^1-0*) = J^ Kfj^,,y{t, s)d ^*^^^^Jt{jt) (the index fi = 



FIG. 3: Dynamic evolution for the initial state d||0) with 
different coupling strength g. The other parameters are T = 
WOmK, ^1 = 2 X lO-'^eV, ^2 = 4 x lO^^el/, ui = uj ^ 
2.5 X IQ-'^el/, CJ2 = 3.5 x 10"^e\/. 



1,2 indicate the left and right baths, and v = b',c' indi- 
cate the fictitious bath b' and c' in the finite temperature 
transformation). Therefore, the exact master equation 
for the double quantum dots system can be derived from 
the SSE, 

dtPr = -i[Hs,Pr] + {Fi{t)[diprj\] + F2 (i) [p,rfi , dj] 



F3{t)[d2Pr,d\] 
F5{t)[dipr,dl] 



F4{t)[prd2,d{] 
Feit)[prdl,dl] 



+ F7{t)[d2Pr,dl] + Fs{t)[prd2,dl] + h.C.}. (18) 

The explicit expressions of the time dependent coeffi- 
cients Fi{t) (i = 1 to 8) can be found in [21|. Here, we 
show some simple properties of this model by plotting 
the time evolution of the population. The detailed study 
of this model will be discussed elsewhere [2J| . In Fig. [31 
we plot the dynamic evolution of the probabilities of all 
the states with different coupling strength between two 
dots. In a long-time limit, the system trends to converge 
to a steady state. When t is small, the electron tunneling 
from one dot to another can be significantly enhanced by 
the direct couplings between the two dots. 

Conclusion. — In this letter, we develop the exact 
stochastic Schrodinger equation approach for solving the 
quantum open system coupled to a fermionic environ- 
ment. The fundamental dynamic equation is derived 
directly from the microscopic quantum model without 
any approximations. By using the Grassmann noise, the 
stochastic Schrodinger equation approach is expanded 
from bosonic to fermionic environments. Three exam- 
ples are presented to show the power of this approach. 
It is worth noting that the stochastic Schrodinger equa- 
tion is versatile enough to deal with a generic fermionic 



environment incorporating cases from a few fermions to 
an infinite number of fermions. The approach can be ap- 
pUed to more reahstic systems when the approximation 
methods are used [24j. 
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